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Abstract 

We consider a simple model for steady-state luminescence of sin- 
gle polymer chains in a dilute solution in the case when excitation 
quenching is due to energy transfer between a donor and an acceptor 
attached to the ends of the chain. We present numerical results for 
Rouse chains without or with hydrodynamic interactions, which are 
taken into account in a perturbative manner. We consider the situ- 
ations of a quiescent solvent as well as the chain in a shear flow and 
discuss the dependence of the steady-state luminescence intensity on 
the strength of hydrodynamic interaction and on the shear rate in the 
flow. 
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1 Introduction 



Luminescent energy transfer in polymers is an important phenomenon. Lu- 
minescent markers are used both for probing the intrinsic polymer dynamics, 
and for probing the properties of the environment using polymers. However, 
the theory of such dynamical phenomena is to no extent satisfactory. The 
problem here is the complicated nonmarkovian dynamics of the system, where 
the most interesting phenomena take place on the time-scales on which the 
systems shows strong memory effects. Even the corresponding initial condi- 
tion problem is hard to solve. No satisfactory quantitative theory exists at 
present for the stationary case. 

Let us start from formulating the problem, and discuss the simplest energy 
transfer model, which will be used throughout the article. Let us assume that 
the ends of a polymer are marked by a donor and an acceptor monomers. 
The molecule is under constant irradiation at a resonant probe frequency, so 
that the donor can get excited with probability A per unit time (we consider 
A as effective intensity of the irradiation). The relaxation of the excited 
state due to spontaneous emission, as well as nonlinear effects connected 
with possible multiple excitation are neglected, so that the only mechanism 
of relaxation is the donor- acceptor energy transfer. We assume that the 
corresponding energy transfer is accompanied by emission of a photon with 
the frequency different from one of the irradiating light. This transfer takes 
place when donor and acceptor approach each other at distance a, hereafter 
called reaction radius. Physically, two situations may take place: Being in 
vicinity of the acceptor the donor still can be excited, and immediately emits 
a photon at the observation frequency. Another situation is the one when, 
being close to the acceptor, the donor gets out of resonance with the probe 
and cannot be excited. In this case, the donor-acceptor system may be in 
one of the two states, on and ojf, being in the on-state the system may be 
excited with probability A per unit time, and emits the photon under the 
transition into the Oj^-state. In what follows the expressions on and off will 
be simply used for denoting states in which the end-to-end distance is above 
and below the reaction radius a, respectively. 

The overall situation might seem simple, however it is much more compli- 
cated than the case of irreversible cyclization [1-7] and is extremely awkward 
for theoretical investigation, even in the absence of flow. Our knowledge 
about the reaction kinetics under flow is sporadic even for simpler reactions, 
see [8,9]. 
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The whole problem would be easily solvable if the life time distributions 
in the on- and the off-siaies were known. Then the probability to be excited 
being in the on-state and therefore the intensity of the emitted light could 
be easily calculated. The probabilities to be in the either state are connected 
with the level-crossing properties of the random process r(i), where r{t) is 
the instantaneous end-to-end distance of the polymer. As for all diffusive 
processes, however, the level-crossing process by r(t) shows a fractal struc- 
ture, so that the mean time between two such crossings is zero (this follows 
immediately from the Rice formula for level crossing density and from the 
form of the two-time correlation function of the end-to-end distances, say, 
for a Rouse polymer, which function lacks the second derivative at zero). 
Again, as for all diffusive processes, this leads to a "tremor" in which r{t) 
crosses the a-level may times until it leaves and performs a long excursion to 
either side. This "tremor" is due to the fact that the diffusion approximation 
(Wiener process) used in the description of the chain (for example through 
the Rouse-like Langevin dynamics) does not adequately mirror a short-time 
dynamics of whatever physical system [10]. However, the existence of this 
theoretical problem does to no extent require for the change of the model 
(say, by introducing underdamped dynamics, as proposed in [10]) since the 
physical problem at hand does not depend on the too-small time behavior 
of the r(t)-process. Indeed, this process is randomly sampled at times ti 
given by a Poissonian flow of photons following with the rate A. The behav- 
ior of r(t) at times much smaller than thus cannot be sampled and can 
physically play no role: this statement is a close analogue of the Nyquist's 
sampling theorem. Thus, the absence of the life-time distributions is not 
a problem of our theoretical model, but a problem of standard mathemat- 
ical tools which rely too much on unphysical, but absolutely unimportant 
short-time properties of a Wiener process. 

Therefore in what follows we concentrate on the numerical investigation 
of the proposed model, and consider the intensity of stationary luminescence 
of the polymer 7(A) under constant irradiation. We discuss the Rouse model 
without hydrodynamic interactions, as well as the role of hydrodynamic in- 
teraction between the monomers, and consider the case when the polymer 
molecule undergoes deformation in a (weak enough) shear flow, which does 
not however cause the full stretching of the molecule. This situation is es- 
pecially interesting as the case when the stationary luminescence of diluted 
polymer solution can be used as a probe for the flow structure. Authors are 
not aware of any experimental realizations of such visualization method, thus 
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our theoretical study might serve as a proof-of-principle for such immediate 
flow diagnostics method. 

2 Simulation approach 

Let us start from discussing our numerical algorithm. Our simulations con- 
sist of two independent parts: the simulation of the r(t)-trajectories, which 
are then stored with high enough resolution, and their analysis giving the 
steady-state luminescence intensity. The reason for this approach is that one 
realization of the process can then be used for getting J(A) for a variety of 
parameters a and A of the model, so that the most time-consuming part of 
the simulation has to be done only once for exactly the time necessary to get 
enough statistics. 

Let us concentrate first of the last part of the problem, namely on the 
evaluation of the stationary luminescence intensity for a given realization of 
r(t)-process. From the record of the r{t) (time resolution of stored data has 
to be much smaller than the minimal A~^ used in simulations) we define 
the a-crossings of the process and, for given a, obtain the lengths of on- 
and Oj(f-intervals, which are ordered and stored. According to the Poisson 
statistics, the probability not to get excited during the on-interval of duration 
ton is exactly exp(— Ato„), thus the probability to emit light after the on- 
excursion is equal to 1 — exp(— Aton)- Since the intensity of emitted light is 
proportional to the overall number of the intervals during which the system 
made a transition into its excited state, we have for the model where the 
off-staXe is not excitable: 

1 n{a,T) 

m = E [1 - exp(-At.)] (1) 

where i numbers the on-intervals, n{a,T) is their overall number, which 
depends on the reaction radius a and on the overall time of simulations T. 
Equation shows that the intervals of very small duration are sampled 
with the probabihty proportional to their lengths so that, as anticipated, 
the fractal structures in vicinity of the concentration points of the level- 
crossings are not resolved and play no role. Using Eq.(IT)) it is possible to 
scan the whole range of intensities A within one run, which is necessary to 
detect nonlinear effects. The situation in which, being in the Oj^-state, the 
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molecule immediately emits light, can be taken into account by adding the 
corresponding intensity to the expression given by Eq.([H). 

/i(A) = /(A) + AP,;^ (2) 

where Pq// is the probability to be in the off-sta.te, i.e. the overall relative 
time spent below a. For example, for the case without flow, it is simply a 
function of relative reaction radius p = a/ ^ {L'^), where {L^) is the mean 
end-to-end squared distance for the chain, 

Poffip) = erf (JIp^ - p^exp (-^p^) . (3) 

Since this simply corresponds to adding a linear function of A to the results 
for the on- off model, we concentrate in what follows only on these results, 
given by Eq.(^. This result holds for all situations without flow. In the 
situation with flow and with hydrodynamic interactions, it is hard to get 
the analytical expression for Poff- The numerical results following from our 
simulations are presented in Tables 3 and 4. 

Let us now turn to simulation of the trajectories. 



2.1 The Rouse model 

We start from the Rouse chain as the simplest model for a polymer [11, 12]. 
A Rouse chain is a set of beads; each one, except for the two end beads, 
is connected to two neighbors by a harmonic potential, so that the overall 
potential energy of the system reads 

N-l 1 

V=Y. i,k\n-n+i\\ (4) 
i=i ^ 

where k is the harmonic spring constant and corresponds to the position 
of the i-th bead. The end beads are connected only to one neighbor. The 
equation of motion of the chain corresponds to overdamped motion under 
the influence of thermal fluctuations: 

ri = — ^ + -Vi, (5) 

7 7 
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where 7 is the friction parameter and rfi is a zero-mean white noise obeying 
the fluctuation-dissipation relation, 



Vnt)vnt')) = '2kBT^^,6^t3S{t - t'). (6) 

In thermal equilibrium, the following relations following immediately from 
the canonical distribution have to hold independently on the model (and are 
always checked numerically as a proof of the quality of the simulation) : 

= l{N-l)kBT (7) 



^ f-tf (8) 

in = ^M^iu^eL, (9) 

where Etot is the total energy, and d and L stand for the bead-to-bead and 
end-to-end distances, respectively. 

We also now apply a shear flow to the system, v = {ay, 0, 0). The shear 
flow is implemented in Eqs.(jHl) by including a term +ayi for the motion in 
the x-coordinate of each bead i, 

^ IdV , ^ , 

ri = — + + a^/i, 0, . 10 

7 OTi 7 

The characteristic intensity of the flow necessary to compare its effects on 
the chain's conformation in different situations is given by the value of the 
dimensionless parameter arR with tr being the Rouse time [7]. 



2.2 Hydrodynamic interactions 

The situation under hydrodynamic interactions is much more involved. The 
standard approaches [13-15] are very accurate but slow, so that we prefer an 
approximate perturbative one. The quality of the corresponding approxima- 
tions is checked by calculating two thermodynamically flxed parameters of 
the chain in quiescent solvent: its mean end-to-end distance and the overall 
energy. We anticipate that especially the end-to-end distance in the chain 
was found to be extremely sensitive to improper incorporation of the hydro- 
dynamic interaction. We conflned ourselves to the situations under which 
the flrst order of the perturbation theory was found sufficient. 
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The hydrodynamic interactions among the beads are modeled within the 
Zimm scheme [16]. The Zimm model is based on the Rouse chain model but 
the equations of motion for different beads are coupled to each other not only 
through elastic forces but also through hydrodynamic forces. Such coupling 
is a long-range one and is introduced through the Oseen tensor [17], that is 
a 3 X 3 tensor defined for each pair of beads (i-j). 



n 



87r?7|f 



Hi 
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(11) 
(12) 



in the direction of fij and 



where X is a unit matrix, f*^- is a unit vector fij/ 

(r'tj^ is its transpose. The viscosity parameter rj can be expressed through 
7 and the bead's size vq since for i — j one has l/GnrirQ = 1/7. Then, 



n 

In what follows we use 7 
reads: 



3rn 



(13) 



47l^iil 

1. The equation of motion for the i-th bead thus 



r' = Zn,[^^n,)- (14) 

The noises ffj acting on different beads are now not independent, otherwise 
the fluctuation-dissipation theorem would be violated. One often writes the 
corresponding equation of motion in the form 

n = Hfi + 2kBTAtf, (15) 

where Ti is the 3A'^ x 3A^ matrix with the diagonal elements being unity (in the 
units where 7 = 1) and with the nondiagonal elements denoting the Oseen 
terms between the corresponding components of velocity of different beads, 
and the matrix A — VH is defined through A • A^ — H. The elements of the 
vector are now independent, zero mean Gaussian white noises. Actually, 
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the computation of the equations of motions in the Euler scheme reads, 



X2{t + At) = X2{t) +Atn fi{t) +^2kBTAtA ip^{t) 



. ^^^^ 

where /,f are the forces due to the harmonic springs for the i-th bead in the 
(3 axis and are the corresponding components of ip. 

The computation of A can be performed exactly by diagonahzing Ti. This 
exact diagonahzation requires an extremely high computational cost for long 
chains. The widely used method based on the orthogonal polynomials de- 
composition (which gives very exact results) is still too slow to get the runs 
long enough for our purposes. Therefore we decided for a simple approx- 
imate approach based on the perturbation expansion of the hydrodynamic 
interaction. 

To do this we write H ets I + rgS, and then expand the square root 
A — y/I + tqS in powers of ro, 



Ap^I + 



+ 



(17) 



Since in the thermal equilibrium the averages {Etot) being the internal en- 
ergy and (L^) (also being a thermodynamical quantity following immediately 
from equipartition) are not modified by the dissipative coupling introduced 
by the Oseen tensor, we can numerically check the validity of the approxima- 
tions for A for different Tq values. We see that {L?) is extremely sensitive to 
incorrect incorporation of the hydrodynamic interaction, and its calculation 
is used as a probe of the quality of the approximation, see Tables 1 and 2. 
The data for Rouse model give us typical error bars for the simulation of the 
exact model on the same scale. 

Looking at the Tables 1 and 2, one can conclude that for ro up to 0.2, 
the second order approximation is sufficient, and for ro up to 0.1, the ffist 
order approximation (much shorter simulations) is accurate enough. In the 
case tq = 0.5 also the second order gets insufficient. Thus, in our simula- 
tions we restrict ourselves to Tq < 0.2. We use a second order Runge-Kutta 
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Table 1: The quality of perturbative approximations for = 21 







< Etot > 


< > 


Theoretical 




30 


20 


Rouse simulation 




30.185 


21.028 




0-order 


29.591 


27.406 


Zimm, ro = 0.1 


1-order 


30.763 


21.020 




2-order 


30.221 


20.181 




0-order 


31.190 


36.923 


Zimm, ro = 0.2 


1-order 


33.114 


23.843 




2-order 


31.676 


21.839 




0-order 


42.058 


73.349 


Zimm, ro = 0.5 


1-order 


49.344 


44.741 




2-order 


53.215 


46.732 



method to solve Eqs. ()16|) with a sufficiently small time step At = 10^^. For 
the results shown in this paper we run 2 ■ 10^ iterations up to a maximum 
time t = 2 - 10^ for a full trajectory needed for adequate statistics. An initial 
thermalization period of 1000 time units is performed in all cases in order 
to start the trajectories from a thermal equilibrium state. The Compaq Al- 
phaServer HPC320 used to run these simulations requires about 3 hours of 
CPU time for iV = 51 when the first order approximation scheme is chosen. 
The second perturbative order requires more than 120 hours of CPU time 
for the same number of iterations and chain length. 

3 Results 

Although the overall role of fiow and hydrodynamical interaction is rather 
clear, the behavior of the intensity as a function of parameters a and ro is 
not trivial. The fiow elongates the molecule, so that the typical end-to-end 
distance grows with a while the hydrodynamical interactions slow-down the 
dynamics of intramolecular relative motion, which increases the characteristic 
time spent in on-state. 

The behavior of /(A) as a function of hydrodynamic radius and fiow 
intensity strongly depends on the relation p between the reaction radius, a, 
and the equilibrium end-to-end distance of the polymer (i.e. the one in the 
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Table 2: The quality of perturbative approximations for = 51 



iV = 51 




< Etot > 


<L^> 


Theoretical 




75 


50 


Rouse simulation 




75.323 


48.122 




0-order 


76.063 


86.674 


Zimm, ro = 0.1 


1-order 


77.412 


54.021 




2-order 


75.926 


47.746 




0-order 


81.399 


130.64 


Zimm, ro = 0.2 


1-order 


83.903 


67.487 




2-order 


79.167 


51.600 



absence of the flow), L = y {L"^)- For p <^ 1 the polymer is typically in the 
on-state, and thus the flow (elongating the chain and making the transition 
into the OjQ^-state less probable) and the hydrodynamic interaction without 
flow (making the change of states slower) work in the same direction and 
lead to the decrease in intensity, as it is clearly seen in FigO 

For p ^ 1 the molecule is typically in the Ojff-state. Increasing flow in- 
creases the probability of switching to the the on-state, and thus leads to 
increase in the steady-state intensity. The hydrodynamic interaction in the 
absence of the flow also leads to increasing the typical time in the correspond- 
ing state. The effects of the flow and the hydrodynamic interactions for p ^ 1 
are depicted in FigO The increasing effect of the hydrodynamic interaction 
has to do with the interplay of two factors. On the one hand, the longer 
on-intervals get even longer under hydrodynamic interaction, and thus give 
smaller contributions to the the overall intensity. On the other hand, increas- 
ing the interaction makes that more shorter on-intervals are now resolved on 
the timescale of A~^, and these contributions in the intensity overweigh the 
loss due to the former of both effects. 

This explanation shows that the role of hydrodynamic interaction is rather 
subtle, and may lead to interesting effects for both regimes (p -C 1 and p ^ 
1), especially when the flow is present. Indeed, the effect of hydrodynamic 
interaction for the cases with a depends in a fine way on all parameters, 
and may act in opposite directions (compare the curves for no flow and high 
flow, ar/j = 6.56, in both panels of Fig^ and the curves for no flow and 
moderate flow, ar^ = 1.05, in the upper panel of Fig|21). 

The values of Poff which are necessary to establish the connection be- 



10 



tween the two situations discussed in the Introduction (Eq.©) are given in 
Tables 3 and 4 for = 21 and for = 51, respectively. The intensities of 
the flows in these tables correspond to the same values of the dimensionless 
flow intensities utr = 0,0.176,1.05 and 6.56 for = 21 and for = 51 
chains (the Rouse times being tr = 13.51 and tr = 84.43, respectively). 

4 Conclusions 

We presented the results of numerical simulations of the intensity of steady- 
state luminescence of single polymer chains in a dilute solution due to excita- 
tion quenching in a simple model in which donor and acceptor are attached 
to the ends of the chain. The chain is modeled by simple Rouse dynamics 
without or with hydrodynamic interactions, which are taken into account in 
a perturbative manner. We consider the situations of a quiescent solvent 
as well as the chain in a shear flow. Depending on the relation between 
the effective distance for energy transfer and the typical end-to-end distance 
of the chain different regimes are encountered with respect to dependence 
of the steady-state luminescence intensity on the strengths of the flow and 
of interaction. Such luminescent probes may be used for experimental flow 
diagnostics. 
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Table 3: Poff for the chain with = 21 



a = 


a = 1 


a = 4 


a = 8 


ro = 


0.014843 


0.502626 


0.974239 


ro = 0.05 


0.014831 


0.499207 


0.972863 


ro = 0.1 


0.013845 


0.484438 


0.969214 


a = 0.0125 


a = 1 


a = 4 


a = 8 


ro = 


0.01467 


0.49977 


0.97302 


ro = 0.05 


0.01481 


0.49775 


0.97191 


ro = 0.1 


0.013724 


0.48283 


0.96835 


a = 0.078 


a = 1 


a = 4 


a = 8 


ro = 


0.013605 


0.443851 


0.937554 


ro = 0.05 


0.013742 


0.456638 


0.947050 


ro = 0.1 


0.012997 


0.450992 


0.949065 


a = 0.488 


a = 1 


a = 4 


a = 8 


ro = 


0.006830 


0.15618 


0.43874 


ro = 0.05 


0.007423 


0.18520 


0.52606 


ro = 0.1 


0.008600 


0.20918 


0.58951 
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Table 4: P^ff for the chain with = 51 



a = 


a = 2 


a = 7 


a = 12 


ro = 


0.034207 


0.613049 


0.960062 


ro = 0.05 


0.033190 


0.594208 


0.957482 


ro = 0.1 


0.028929 


0.556314 


0.945990 


a = 0.002 


a = 2 


a = 7 


a = 12 


ro = 


0.03386 


0.61418 


0.95983 


ro = 0.05 


0.03270 


0.59412 


0.95756 


ro = 0.1 


0.02932 


0.55641 


0.94633 


a = 0.0125 


a = 2 


a = 7 


a = 12 


ro = 


0.033417 


0.575235 


0.933795 


ro = 0.05 


0.030730 


0.571878 


0.945019 


ro = 0.1 


0.027587 


0.540211 


0.939310 


a = 0.078 


a = 2 


a = 7 


a = 12 


ro = 


0.01609 


0.20126 


0.45965 


ro = 0.05 


0.02417 


0.28267 


0.61573 


ro = 0.1 


0.01273 


0.30860 


0.67650 
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Figure 1: The intensity of steady-state luminescence in the on-Oj^ model as 
a function of irradiation intensity A for a reaction radius such that p <^ 1. In 
panel (a) N — 21 and a — 1. In panel (b) = 51 and a — 2. In both panels 
the same notation is used: the symbol indicates the value of the bead radius 
Tq: filled squares (0). empty circles (0.05) and starts (0.1), whereas solid, 
dotted and dashed lines correspond to aTR = 0, 1.05 and 6.56, respectively. 
tr — 13.51 for the chain with 21 beads and tr — 84.43 for chains with 51 
beads. 
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Figure 2: Same as in Fig^ but now for p ^ 1. In panel (a) TV = 21 and 
a = 8. In panel (b) N = 51 and a = 12. We use the same notation for the 
lines as in FigQ 
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